Constraints on Neutrino Mass from Sunyaev— Zeldovich Cluster Surveys 
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The presence of massive neutrinos has a characteristic impact on the growth of large scale struc- 
tures such as galaxy clusters. We forecast on the capability of the number count and power spectrum 
measured from the ongoing and future Sunyaev- Zeldovich (SZ) cluster surveys, combined with cos- 
mic microwave background (CMB) observation to constrain the total neutrino mass M,^ in a flat 
ACDM cosmology. We adopt self-calibration for the mass-observable scaling relation, and evaluate 
constraints for the South Pole Telescope normal and with polarization (SPT, SPTPol), Planck, and 
Atacama Cosmology Telescope Polarization (ACTPol) surveys. We find that a sample of ~ 1000 
clusters obtained from the Planck cluster survey plus extra information from CMB lensing extrac- 
tion could tighten the current upper bound on the sum of neutrino masses to = 0.17 eV at 68% 
C.L. Our analysis shows that cluster number counts and power spectrum provide complementary 
constraints and as a result they help reducing the error bars on by a factor of 4 — 8 when both 
probes are combined. We also show that the main strength of cluster measurements in constrain- 
ing My is when good control of cluster systematics is available. When applying a weak prior on 
the mass-observable relations, which can be at reach in the upcoming cluster surveys, we obtain 
c"A/„ = 0.48 eV using cluster only probes and, more interestingly, am^ = 0.08 eV using cluster -|- 
CMB which corresponds to a S/N ~ 4 detection for > 0.3 eV. We analyze and discuss the de- 
generacies of My with other parameters and investigate the sensitivity of neutrino mass constraints 
with various surveys specifications. 



I. INTRODUCTION 

Measuring masses of neutrinos is one the major goals 
of particle physics and cosmology. While atmospheric 
and solar neutrino oscillation experiments are sensitive 
to neutrino flavor, mixing angle, and the mass difference 
among different species, cosmological data are instead 
more sensitive to the absolute mass scale M^, = ^ m^. In 
fact, the most stringent upper bound of the total neutrino 
mass is coming from CMB and large scale structures since 
massive neutrinos leave detectable imprints throughout 
the history of the universe. Most recently, [l| obtained 
< 0.44 eV at 95% C.L. by combining CMB from 
WMAP 9 year data, BAG from Sloan Digital Sky Survey 
(SDSS) and WiggleZ, and Hq data from Hubble Space 
Telescope (HST). 

In this work, we explore the prospects of employing on- 
going and future galaxy cluster surveys detected by the 
SZ effect in constraining neutrino mass. Galaxy clus- 
ters are in principle a powerful tool for probing neu- 
trino properties. Neutrino becomes nonrclativistic af- 
ter the epoch of decoupling if its mass scale is smaller 
than 0(0.1) eV. The relativistic behavior of neutrinos, 
as opposed to cold dark matter, causes the suppression 
of matter perturbations on small scales with respect to 
the case in which neutrinos are massless and all dark 
matter is cold. The presence of massive neutrinos affects 
the growth rate of perturbations in the linear regime, 
and, as a consequence, the shape of the matter power 
spectrum and cluster abundance. Current measurements 
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from X-ray cluster surveys obtained a tight upper limit of 
Ml, < 0.33 eV by combining measurements of Chan- 
dra X-ray observations of galaxy clusters, CMB from 
WMAP 5 year data, BAG and type lA supernova (both 
from HST Key Project). Similarly, measurement from 
galaxy power spectrum from the SDSS-III BAG survey 
+ CMB + SN found < 0.34 eV 0. 

Subsequently, several works were dedicated to discuss 
the prospects of utilizing large future surveys of large 
scale structures (galaxy or galaxy clusters) in different 
wavelengths (e.g. |4|-|ll|). These works showed that con- 
straints of neutrino mass depend on assumptions of the 
underlying cosmology (e.g. inclusion of dark energy or 
flatness), cluster physics, and the use of external priors 
(e.g. CMB lensing extraction) . Here we revisit the analy- 
sis to forecast the constraint of the total neutrino mass, in 
the framework of flat ACDM universe and, like past con- 
straints, the standard scenario with only three neutrino 
species. We use cluster abundance and power spectrum 
as the observables that will be obtained from various SZ 
cluster surveys: the Planck, ACTPol, SPT, and SPTPol 
cluster surveys. These surveys are very promising and, 
in the next couple of years, will provide large samples 
of mass selected clusters out to high redshift. With re- 
spect to previous works [e.gH, [l3| which also employ SZ 
cluster surveys, we provide a more realistic survey specifi- 
cations to characterize the cluster detection and include 
the self-calibration to characterize the uncertainties of 
the mass-observable relations. We also discuss the de- 
generacy of the neutrino mass with dark energy, which 
is lacking in previous studies, and compare the strength 
of cluster probes with CMB on the constraining power of 
neutrino mass. 

The paper is organized as follows. In Sec.|ll]we discuss 
the effects of M^, on the large scale structures. In Sec. IIIII 
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we present the methodology which includes the descrip- 
tion of the future SZ cluster samples and the fisher matrix 
formalism. The main results are presented Sec. IIVI and 
discussed in Sec. El Finally, a conclusion is presented 
in Sec. ED 



II. IMPACT OF NEUTRINO MASSES ON 
GROWTH OF THE LARGE SCALE 
STRUCTURES 

The presence of massive neutrinos mildly affects ex- 
pansion history but significantly impacts the growth of 
structure through free-streaming. Fluctuations on co- 
moving scales that enter the horizon when neutrinos are 
still relativistic may be reduced in amplitude because 
neutrinos would tend to leave the perturbation. This 
effect, that is neutrino mass dependent, typically oc- 
curs on scales below the free-streaming scale: kfg ~ 
l.b^yriMh^/il + z){M^ /eY)M-pc-\ Thus, the growth of 
any structure that have wavenumber larger than fcfg will 
be less efficient. A smaller neutrino mass increases the 
free-streaming scale, but also reduces the neutrino frac- 
tion with respect tot he total amount of dark matter, 
mitigating the overall suppression. 

As a result of these dependences measurements of the 
large scale structures such as cluster number counts and 
power spectrum can be used to place constraints on neu- 
trino masses. 

The late-time evolution of perturbations in a ACDM 
cosmology with massive neutrinos can be accurately de- 
scribed by the product of a scale dependent growth func- 
tion and a time dependent transfer function. For ex- 
ample, [l2j derived a reasonable approximation to the 
analytical expression of the transfer function for small 
scales. In this work, we employ the transfer function de- 
termined numerically from CAMB [l3| which provides 
precise estimate on the matter power spectrum and in- 
clude non-linear effects at large-fc limit in which the an- 
alytical expressions fail to give an accurate estimate. 

III. ANALYSIS 

Our analysis closely follows the treatment of , here 
wc only outline the method and refer the readers to 
for details. For cluster abundance and clustering, we use 
the results of numerical simulations from [l5| for the clus- 
ter mass function n{M, z), and [l^ for the halo bias. 

A. Cluster survey 

We consider five upcoming SZ cluster surveys studied 
in : the Planck survey, the South Pole Telescope nor- 
mal and polarization survey (SPT and SPTPol respec- 
tively), the Atacama Cosmology Telescope polarization 
survey (ACTPol). Each of these surveys has different 



TABLE I: Properties of SZ cluster survey 



Survey 


Area (sq. deg) No. of clusters 


Planck 


30000 


2500 


SPT 


2500 


500 


SPTPol 


625 


1000 


ACTPol 


4000 


500 



specifications for the selection threshold, i.e. Miim(z), 
and their properties are summarized in Tab. ID 

Briefly, for the Planck survey we adopt a flux limit 
of i200,pc > 2 X 10~^arcmin^ where >200,pc is the 
integrated comptonization parameter within the radius 
enclosing a mean density of 200 times the critical density. 
This corresponds to a 5(t detection threshold and would 
yield ~ 1000 clusters. 

For the SPT survey, i.e. single frequency at 150 GHz, 
wc employ the calibrated selection function of the survey 
by [i^ and adopt a detection threshold at 5a. This yields 
~ 500 clusters . The SPTPol has an increased sensitivity 
at 150 GHz than the normal survey and we account for 
this, following previous work, by scaling the mass limits 
by a factor of 3.01/5.95. The expected number of clusters 
is - 1000. 

For the ACTPol survey, we include clusters with 
^-^^200,^-0 > 5 X IO^^Mq/i"^ (Sehgal 2011, private com- 
munication) which corresponds to a 90% completeness. 
This straight mass cut resultsq in ~ 500 clusters. 

We construct cluster sample for the Planck survey in 
the redshift range < z < 1. We impose a lower cut 
Zeut = 0.15 for the SPT, SPTPol, and ACTPol survey. 
Currently, the SPT team is setting a low redshift cut at 
Zcut — 0.3 in their released cluster sample, due to difficul- 
ties in reliably distinguishing low-redshift clusters from 
CMB fluctuations in single frequency observations. Nev- 
ertheless, with upcoming multi-frequency observations, a 
lower cut Zcut = 0.15 will likely be attained. We therefore 
apply this cut in our work. 



B. Fisher matrix forecasts and cosmological 
parameters 

We estimate the constraints on cosmological parame- 
ters by applying the Fisher matrix formalism to future 
SZ cluster surveys. This approach can best approximates 
the likelihood when the fiducial model is close to the true, 
underlying model and the likelihood is close to gaussian. 
Typically, the gaussian approximation is more accurate, 
and the use of the fisher matrix better justified, when the 
likelihood is peaked and the parameter in hand has little 
degeneracies with other parameters. In order to achieve 
this goal, the use of external priors an be beneficial. 

For example, [l^ noted that the CMB power spec- 
tra likelihood function for the neutrino mass differs from 
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the gaussian case due to strong parameter degeneracies, 
particularly for models with many parameters. These 
authors suggests the use of CMB lensing extraction in- 
formation in order to sharpen the likelihood and make 
it better approximated by Gaussian. We will adopt the 
same strategy, as described below. 

The Fisher matrix for the cluster number counts and 
power spectrum is described in detail in [T^]. Similarly, 
for our main results, we here consider self-calibration to 
account for the uncertainties of the observed cluster mass. 
We add the Planck CMB lensing extraction (LE) that 
is considered to be a very promising way to constrain 
neutrino mass (e.g. [HI). The CMB anisotropics obey 
Gaussian statistics in the absence of weak lensing, and 
therefore they are fully described by the temperature and 
polarization power spectrum. Weak lensing, however, in- 
troduces non-gaussianity in both the temperature and 
polarization anisotropies [lO, [llj. Therefore, extracting 
the lensing information from CMB (e.g. using quadratic 
estimators [2^ - l2H ) would provide the lensing potential 
and delensed CMB anisotropies, and hence extra infor- 
mation to the fisher matrix. In the following, we refer 
to the fisher matrix results obtained from CMB lensing 
extraction as the CMB LE. As shown in [11], CMB LE is 
useful in providing strong neutrino mass constraints and 
potentially breaking of the major neutrino mass degen- 
eracies with other parameters [l^. While very promis- 
ing, the exploitation of higher order statistics may suf- 
fer from subtle ways from the effect of galactic and ex- 
tragalactic contaminants. For this reason, we also con- 
sider constraints coming from the CMB power spectrum 
only (with lensing) when combining probes with cluster' 
s ones. 

We adopt a spatially flat ACDM model as the fiducial 
model. The set of parameters included in our analy- 
sis is (r2f,/i^, r2j\//i^, J7a, Mi^, rig, (787 w^o, "u^a)- The fiducial 
values arc adopted from the best fit flat ACDM model 
from WMAP 7yr data, BAG and Hq measurements [l^] : 
rith'^ = 0.0245, riMh'^ = 0.143, = 1 - fin = 0.73, 
AU = 0.3 eV, Us ^ 0.963, as = 0.809, wq^ -l,Wa^O. 

As proposed in [H, [2^, we can use cluster surveys 
to constrain the mass observable relation by considering 
self-calibration, hence taking into account the systematic 
errors of the SZ surveys due to uncertainties in observed 
cluster mass. In this work, we follow [30j to introduce 
four nuisance parameters, Bmo^ Q!, crinA/.Oi that spec- 
ify the magnitude and rcdshift dependence of the frac- 
tional mass bias Bf,[{z) ~ i?7\/.o(l + z)°' and the intrinsic 
scatter crinA/(-z) = crinA/,o(l + 2)^. We adopt fiducial val- 
ues of i?A/o = 0, a = 0, CTinAf,o = 0.1, /3 = 0, hence 
corresponding to zero mass bias and 10% intrinsic scat- 
ter. In deriving the main results, we will not make any 
assumption on the four nuisance parameters and leave 
them free to vary. We discuss the impact of this assump- 
tion in Sec. HVCl 



IV. RESULTS 

A. Cluster number count and power spectrum 

Tab. m] summarizes the neutrino mass constraints from 
the Fisher matrix analysis for Planck CMB (with and 
without LE), cluster number counts, and power spectrum 
for the four cluster surveys. Constraints of My from clus- 
ter number counts alone are better than power spectrum 
ones, however, each of them is very weak when considered 
separately, with (Tm„ > 4 eV in the context of ACDM 
model. When combining information from both probes, 
the constraints are improved significantly by a factor of 
4 — 8. The best case is obtained from the Planck cluster 
survey with fT7\/_^ = 0.94 eV, whereas the constraints from 
other surveys are a factor of two worse. 

B. Cluster probes + CMB 

Adding the Planck CMB priors breaks degeneracies 
(see Sec. IV A[) and improves the constraints (number 
count or power spectrum alone) further by a factor of 
> 4 (without LE) and > 5 (with LE). When including 
all the information but LE, i.e. count -I- power spec- 
trum + CMB, we find the best constraint comes from 
the Planck and ACTPol cluster survey with aM„ =0.23 
eV. This is 80% better than that obtained from Planck 
CMB alone (cta/^ 0.41 eV). Including CMB priors also 
shrinks the difference in gm^ among different surveys in 
which it is now = 0.23 — 0.30 eV. Similar results are 
obtained when we add the CMB LE and the best con- 
straint is ctm^ = 0.17 eV. This suggests that the improve- 
ments in fjM^ are mainly driven by CMB information. 
We note that a perfect cleaning of all the astrophysi- 
cal foregrounds is assumed when computing the CMB 
fisher matrix in this work. Foreground contamination 
dominates at small angular scales (e.g. / > 1000) and 
would introduce extra non-gaussianity and spoil the lens- 
ing extraction process [sij. Nevertheless, [Iq found that 
the effect of no foreground subtraction in Planck CMB 
(with and without LE) only degrades the constraint 
marginally (by 9%). Therefore, our results that involve 
CMB information can be considered to be robust against 
foreground contamination. 

We repeat the analysis with a fiducial Af^ = 0.1 eV 
instead to investigate the effect on the constraint with 
less massive neutrinos. The results are very close (within 
15%) to those for = 0.3 cV when using cluster probes 
only, and are almost the identical when CMB priors are 
added. 



C. Self-calibration and uncertainty of nuisance 
parameters 

The dominant systematic errors for SZ derived con- 
straints are the uncertainties in the mass observable re- 
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TABLE II: Marginalized la errors on Af^ (in units of eV). 



CMB CMB LE 


Survey 


no prior 
NC P{k) Comb 


+ CMB prior 
NC P{k) Comb 


+ CMB LE prior 
NC P{k) Comb 


fiducial = 0.3 eV 


0.41 0.21 


Planck 
ACTpol 
SPT 
SPTpol 


4.06 7.83 0.94 
6.04 11.73 3.33 
12.44 12.45 2.12 
12.59 7.81 1.79 


0.29 0.29 0.23 
0.29 0.29 0.23 
0.33 0.31 0.30 
0.32 0.30 0.28 


0.20 0.20 0.17 
0.20 0.20 0.17 
0.21 0.20 0.20 
0.21 0.20 0.19 


fiducial = O.f eV 


0.52 0.19 


Planck 
ACTpol 
SPT 
SPTpol 


5.09 19.98 0.77 
17.97 48.63 2.83 
6.56 31.39 1.78 
12.59 7.81 1.79 


0.37 0.42 0.23 
0.43 0.38 0.26 
0.43 0.39 0.28 
0.32 0.30 0.28 


0.18 0.19 0.15 
0.19 0.18 0.16 
0.19 0.18 0.17 
0.21 0.20 0.19 



lation due to structure and evolution of clusters. We 
can ask how much could be gained by eliminating such 
uncertainties. For example, we can expect some exter- 
nal constraints on the nuisance parameters by using de- 
tailed studies of individual clusters or combining differ- 
ent information from optical, weak lensing. X-ray and SZ 
measurements. To estimate the effect of self-calibration 
of systematic uncertainties on the neutrino mass con- 
straints, we repeat the forecasts with different priors on 
the four nuisance parameters as summarized in Tab. IIIII 

We first discuss the results when applying a "weak" 
prior, i.e. using current knowledge on the calibration 
on the mass proxies Acta/.o = 0.1, A/3 ~ 1, A_Bm,o ~ 
0.05, Aa = 1. In the case of cluster count + power 
spectrum, the la error reduces marginally for SPT and 
SPTPol, but significantly (by a factor of two) for Planck 
and ACTPol. This results in ctm^ = 0.48 eV for the 
Planck cluster survey which is competitive with the CMB 
only constraint. In the case of adding the CMB (with and 
without LE) priors, the Itr errors generally reduce by a 
factor of two and resulted in, for the best case as obtained 
by the Planck survey, (Tm„ = 0.08 eV, which corresponds 
to a S/N « 4 detection for > 0.3 eV 

Similar results are obtained when applying a "strong 
prior" , i.e. the four nuisance parameters are held fixed 
at their fiducial values, which is equivalent to assuming a 
perfect knowledge of cluster true masses. The constraints 
are improved significantly by 66 — 236% in the case of 
cluster count -I- power spectrum, and a factor 2 — 3 when 
the CMB priors are further added. The best constraint 
is, again with the Planck cluster survey, cfm^ = 0.07 eV 
which is a relative marginal improvement with respect 
to the weak prior case. While it is unrealistic to have 
perfect knowledge on the mass observable relations, one 
can achieve similar scenario by restricting the analysis 
to a relatively small subset of clusters for which follow 
up observations are available. This would ensure a sam- 
ple with well calibrated mass proxies. For example, it 
has been shown in Q that the ability to constrain dark 
energy parameters from a small sample of « 50 well cal- 



ibrated X-ray clusters is comparable to a larger sample 
of « 10000 optical clusters (e.g. SDSS 0). 

Unlike other parameter constraints (e.g. non- 
Gaussianity with galaxy clusters [l^l ) , the results of the 
weak prior are sufficiently close to the those from the 
strong prior. The prospect of achieving the weak prior 
conditions is promising, e.g. clusters detected in weak 
lensing measurements or a subsample of objects hav- 
ing extensive multi- wavelength follow-up. Therefore the 
cluster probes are good enough to provide interesting 
constraint even without perfect knowledge of the scaling 
relations. 

As a final remark, we would like to compare our 
count -I- CMB result with [l^ which similarly presented 
constraints assuming perfect knowledge of cluster 
mass and used Planck cluster count -I- CMB. Our result 
(fM„ = 0.17 eV) is a factor of 2.8 worse than that ob- 
tained in [l^ . We note that the discrepancy is due to the 
different assumption on the total number counts: ~ 6000 
in (loj and « 1000 in this work for the Planck survey if 
a 5f7 survey detection limit is assumed. Our estimate is 
based on the conservative assumption that ensures high 
level of completeness (90%) and realistic mass limits that 
vary at different redshifts, while [l3| assumed a constant 
and lower mass threshold. 



V. DISCUSSION 

A. Parameter Degeneracies 

The dark energy equation of state wq and Mi, is one of 
the major parameter degeneracies. Fig. [T] shows the la 
constraints on and Wq computed from cluster number 
counts, power spectrum, combination of the two, with 
and without LE of the Planck CMB. The contour for 
number count shows a clear diagonal alignment, and the 
degeneracy direction can be understood as follows: an 
increase in neutrino mass suppresses the growth of struc- 
ture formation, this can be compensated by a larger rate 
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TABLE III: Fractional improvement "i^iu.no with various priors (see Sec. IIV CI) . Tire values of am,, no of the corre- 

'^M„,uj<!ak/ strong ^ ' 

sponding cases are those from Tab. [Til A fiducial — 0.3 eV is assumed. The best case is obtained by the Planck survey, 
aA/„ = 0.08 eV (same for + CMB or + CMB LE). 



Probes 


prior 


Planck ACTpol SPT SPTpol 


dN/dz + P{k) 


weak 


1.97 


2.43 


1.18 1.06 




strong 


2.34 


3.36 


1.67 1.66 


dN/dz + P{k) + CMB 


weak 


2.76 


2.16 


2.12 2.53 




strong 


3.25 


2.53 


2.33 2.88 


dN/dz + P{k) + CMB LE weak 


2.25 


1.67 


1.64 1.95 




strong 


2.63 


1.92 


1.77 2.19 



of accelerated expansion (i.e. more negative w). The 
constraints from power spectrum is less degenerate but 
show different degeneracy directions. As a result, com- 
bining information from both probes greatly improve the 
constraints. To sec the effect of Wq on constraint, we 
derive crj\/_^ again by marginalizing over wq and Wa- Wc 
find that, as expected, only the constraints from number 
count are affected (improve by a factor of > 2), while 
those from power spectrum are barely affected. Further- 
more, only modest improvements are obtained when com- 
bining number count and power spectrum in this case. 

The degeneracy between curvature fi^ and neutrino 
mass is also known to be significant and impact on 
both Mu and the number of neutrino species N^f / , which 
could affect the constraints coming from CMB [e.g [H, 
Is^ l . However wc note that the cluster probes used in this 
work are related to the growth of structures which are not 
sensitive to ^k- Thus we expect that including ^Ik in 
the fisher matrix analysis would not impact our results. 
It is out of the scope of this paper to study in depth 
the impacts of including an extended set of parameters 
(e.g. r^K, Neff). Nevertheless it would be potentially 
interesting to study their effects for growth of structures 
and we leave it for future works. 



B. Survey sensitivities to neutrino mass constraints 

In order to better understand what aspects of SZ sur- 
veys would improve the constraints on neutrino mass, wc 
repeat the fisher matrix calculation that includes differ- 
ent range of wavenumber k and cluster mass M. 

One of the major dependences is the maximum k val- 
ues (fcinax) as it determines the smallest scales that can 
be probed by a survey. The effect of massive neutrino is 
particularly prominent at small scales (large k values), in 
which free streaming of neutrinos prevent structure for- 
mation. In this work, we use the same fcmax = 0.1 h/Mpc 
in the power spectrum fisher matrix for all surveys con- 
sidered. We do not attempt to increase the fcmax beyond 
this value to avoid the non-linear effects at smaller scales. 
Furthermore, this scale is the limit that can be reached 
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FIG. 1: Joint constraints on the Mi, and wq. All curves de- 
note 68% confidence level, and are for number counts only 
(blue), power spectrum only (cyan), and combination of the 
two (green), Planck CMB (yellow), and Planck CMB LE (dot- 
ted yellow). 



by SZ cluster surveys. Instead we study the dependence 
when small scale modes are lost, as shown in Fig. [2] (left). 
The effect of losing small scales information begins at 
k « 0.06 h/Mpc, which corresponds to the free stream- 
ing scale fcfs at z = 2, with Mi, — 0.3 eV. The prospect of 
using smaller scale modes in constraining neutrino mass 
would be coming from galaxy surveys which can probed 
down to fc w 0.5 Mpc/h. A number of studies forecasted 
the neutrino mass constraints from future galaxy sur- 
veys (e.g. [13 (BOSS) and (EUCLID-like)) and found 
that the improvement in (Tm„ beyond fc = 0.1 Mpc/h is 
marginal (sec Fig. 6 in HH). 

The other relevant dependence is the limiting cluster 
mass that differentiates the various SZ surveys. We show 
in Fig. [5] (right) the fully marginalized (Tm„ from the 
power spectrum -I- Planck CMB as a function of min- 
imum cluster mass used in the calculation. It is clear 
that a deeper survey that can probe down to lower mass 
would improve the constraints. It is also interesting to 
note that the SPTPol survey, despite of its small sky cov- 
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erage, can perform better than the ACTPol and SPT sur- 
vey because it can detect clusters down to « 2 x 10^"* Mq. 
Therefore a deep survey can compensate for the sky cov- 
erage when constraining neutrino mass. 

VI. CONCLUSION 

In this work, we explored the possibility of using future 
and upcoming SZ cluster surveys to constrain neutrino 
masses. We employ the Fisher Matrix analysis to forecast 
the sensitivities of various SZ surveys in constraining the 
total neutrino mass Mi, in the context of flat ACDM 
cosmology. We do so by making use of the cluster number 
counts and power spectrum, and taking into account the 
self-calibration of mass-observable scaling relations. 

In general, we find that the M^, constraints from clus- 
ter number count and power spectrum is weak if they 
are considered separately, due mainly to strong parame- 
ter degeneracy between M^, and wq, especially in the case 
of cluster number count. However, such degeneracy can 
be broken if the two probes are combined, which helps 
to improve the constraints considerably. For example, a 
sample of ~ 1000 clusters obtained from the Planck clus- 
ter survey gives ctm^ =0.94 cV (0.43 eV with weak prior). 
The constraints can be further improved when combined 
with CMB priors. The best constraint is obtained for the 
Planck and ACTPol survey, with aM„ = 0.23 eV (CMB) 
and crA/„ = 0.17 eV (CMB LE). This is « 80(25)% im- 
provement with respect to the CMB (CMB LE) only con- 
straint. The use of CMB Icnsing extraction can better 
help the cluster only constraints because it determines 
the neutrino's free streaming effect on the matter power 
spectrum and break some of the parameter degeneracies. 

While we find that Al^, constraint is mainly driven by 



CMB and the addition of cluster probes, i.e. number 
count -I- power spectrum, to CMB only helps marginally, 
the use of clusters is still beneficial if we have good con- 
trol of cluster systematics. For example, when applying a 
weak prior on the mass-observable relation, the la error 
on Ml, as obtained from cluster count + power spectrum 
goes down to 0.48 eV and is competitive with CMB only 
constraint. If wc further combine with CMB priors, 
reduces to 0.07 eV, which corresponds to a w 4cr de- 
tection for Ml, > 0.3 eV. The prospect of achieving the 
weak prior conditions is promising, e.g. clusters detected 
in weak lensing measurements or a subsample of objects 
having extensive multi-wavelength follow-up. Therefore, 
cluster measurements are useful, as an independent probe 
of the Mij with respect to the CMB, in tightening the 
current bound on Mi,. 

Wc find that a deeper cluster survey that detects 
smaller mass clusters, e.g. down to 2 x IO^^A/q like 
SPTPol, improves neutrino mass constraints. This is 
because of the effect of free streaming of massive neu- 
trinos that prevents structure formation to happen at 
small scales. Likewise, the availability of the small scale 
modes, i.e. the maximum k values that can be probed 
by a cluster survey, also helps the constraints. We show 
that the modes at fc > 0.06 h/Mpc are important as they 
help decreasing ctm^ significantly. 
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